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a physical system. 

A UP |UOT*»A1* luivuuw** *w*^»w« R " 

Such processing may be used m data Aversion and a particular example of this is the 
5 inversion of time-lapse seismic data. 

Seismic reflection is a technique used to determine details of structures beneath the 
surface of the Earth. The resolution that may bo achieved makes mis technique the 
method of choice fox oil exploration and mapping of subsurface rock structures. It is 
10 also applicable to experimental research that probes the fine structure within the Earth's 
crust and at the crust-mantle boundary. 

the technique involves generating downward-propagating seismic waves in succession 
at a number of locations, within the region being explored. A large number of receivers 
are positioned at intervals away from each source location and these receivers record the 
amplitudes (for example, in terms of pressure, displacement or its derivative) of senmuc 
waves reflected back up to the surface from subsurface inhomogeneities over a period 
of time. The recorded waves are usually deconvolved, removing the effects of the 
source and receiver (which have their own response functions). 

Reflection data typically have low arrrplitudes and are contaminated by multiple 
reflections and other kinds of noise. Various acquisition and processing techniques may 
be used to improve dgnal-to-noiae ratios, such as averaging (stacking) of traces w,m the 
Bame midpoint, taking into account different distances between source and receiver, and 
25 discrhmnation of multiple reflections based on either their periodicity or wavefkont 
angles which differ from the primary reflections. Further, the data may be cc™* 
positioned in space by a process called migration, which moves dipping events into 
their correct position. When comparisons are made between two or more datascts over 
the same area, careful analysis between the amplitude, time and other attributes of the 
30 datasete may be made. 



20 
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After the appropriate crarections- which may fiirther include correction for other known 
ratal variables, the data axe combined to provide a graphical representation of 




u* in* 



die subsurface i 



tliTtinv* 



geneities. 



5 Seismic reflection data obtained by field experiments are then processed to obtain a 
three dimensional image of subsurface structure as described above. The three 
dimensions refer to the spatial dimensions "illuminated" by the seismic data. The 
vertical axis may represent depth, or two-way vertical seismic wave travel time. 

10 The amplitudes of reflected seismic waves are indicative of the subsurfeee mflection 
strengths, contaminated by noise. The reflection strength depends upon the reflection 
coefficient, which may be defined as a function of the relative contrasts of the elastic 

of flie subsurface layers. 




15 The clastic properties of an isotropic, elastic medium are completely described by three 
parameters, for example the two elastic Lame parameters and the density. Other 
panmeterisations are possible, for example acoustic impedance, shear impedance and 
density. A third example is P-wave velocity, S-wave velocity, and density. The 
transfonnation between different sets of elastic parameters is well defined and 

20 straightforward- 

In general, me elastic properties vary spatially. In order to explain the relationship 
between the elastic properties and the seismic data it may be convenient to imagine the 
subsurface as a stack of geological layers. The layer properties are described by the 
25 elastic properties of the rocks within the layers while the seismic data are related to the 
contrasts of the layer properties between successive layers. The seismic data are 
therefore suitable for interpreting subsurface layer structures since they image the 
boundaries between the layers. 

30 Seismic inversion is defined herein as the process of t_forming (inverting) seismic 
reflection data to elastic material properties, i.e. taking amplitudes (measurements of 
contrasts) and using them to infer physical layer properties. Numerous mflerent seismic 
inversion techniques are known. • 

♦ 
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Over a period of time, certain types of rock, known as source rocks, will produce 
hydrocarbons. The produced hydrocarbons are then transferred to and stored in rocks 
known as reservoir rocks through various geological processes. During production of 
hydrocarbons in a subsurface region, the effective elastic material properties of the 
reservoir rocks change with production time, where production time is the fourth 
dimension in seismic 4D analysis. The changes of the effective elastic properties of the 
reservoir rocks may be caused by changes of the pore fluid saturations in the reservoir 
rooks, but also by pressure and temperature changes. Explained by a simple layer-based 
earth model concept, the properties of the reservoir layer are changed during production, 
implying changes in me reflectivity for the upper and lower reservoir interfaces. The 
measurements taken at a further seismic survey are related to me new contrasts at the 
boundaries between adjacent layers. 



IS 



20 



25 
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Reservoir changes are often interred nom a w>"P«« v " — — > ~ 

amplitudes of seismic waves reflected at interfaces bounding or within the reservoir) for 
Afferent seismic surveys acquired at different stages of the production. A more direct 
interpretation can be based on difference data. Difference data are established by 
subtracting two time-separated seismic surveys covering a common part of the earth. 
The difference data, after me proper time-alignment during pre-processing represent a 
spatial image of me changes of the relative contrasts between the two different 
acquisition times. 



For a three dimensional seismic dataset, the classic inversion problem is to estimate the 
elastic material parameters from the three dimensional seismic data. A natural 
extension of 3D inversion to inversion of time-lapse seismic data (4D) is to invert &e 

different 3D dataseto separately by a known method, and then subtract the results to 
obtain the changes* 

However, me reliability of 4D interpretations is difficult to assess, and are made by 
qualitative assessment. A full consideration of me uncertainties involved is important 
for making an accurate inference of the changes in the reservoir properties between the 
two seismic surveys. The results of such seismic analysis maybe important in reservoir 
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management in that the inferred reservoir properties are used to evaluate, for example* 
new drilling targets and future drainage strategies. 



10 



15 



25 



30 



inversion provides quantitative estimates of the elastic reservoir properties. 
However, inversion of noisy seismic data is known to be a difficult and ill-posed 
procedure. An appropriate assessment of the uncertainties in 4D inversion data has not 
previously been possible. 

Commercial time-lapse inversion techniques have become available, but only with brief 
descriptions of the method Sow? results have been published (Mesdag et al, 2003, 
Integrated AVO reservoir characterisation and time-lapse analysis of the Widuri field, 
65 th Mtgn Bur., Assn. Expl. Geophys., Extended Abstracts). Such methods apply 
separate inversions of the data with some constraint between the results, e,g. a common 
background model. The time-lapse change is then calculated from the change in 
inverted parameters, Sarkar et al, 2003, On the Inversion of time-lapse seismic data, 
73 rd Aim. bitemat Mtg.: Soc. Of Expl. Geophys., U8JM492> mentions inversion of 
seismic differences, but provides no detail of tttf implementation. None of these 
inversion techniques provide uncertainty bounds on the results. 



20 According to a first aspect of the invention, there is provided a method as defined hi the 
appended claim 1. 

Further aspects and embodiments of the invention are defined in the other appended 



claims. 



It is thus possible to provide a technique which permits improved inversion of seismic 
data representing a physical system. Such a technique may be used to handle errors 
intrinsio to such data and can provide, for example, probability distributions or 
uncertainty bounds on the results of inversion. 

For a better understanding of the present invention and in order to show how the same 
may be carried into effect, preferred embodiments of the invention will now be 
described, by way of example, with reference to the accompanying drawings in which: 
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Figuio 1 is a flow diagram illustrating si method constituting an embodiment of the 
present invention; 

Figures 2a to 2c illustrate a reservoir model; 
5 Figures 3a to 3d illustrate modelled seismic data fox first and second surveys and their 

associated noise characteristics; 

Figures 4a to 4c and 5a to 5c illustrate the results of an example of the method of Figure 

1; 

Figure 6 illustrates the uncertainties in the results illustrated in Figures 4a to 4c and 5a 
10 to 5c; 

Figure 7 is a probability map derived ftom the results of Figures 5a to 5c; and 

Figuw 8 ia a block schematic diagram of an apparatus for performing the method of 

Figure 1. 

The embodiments of the present invention relate to a method of data inversion that 
operates directly on seismic difference data, and in particular to the difference between 
two sets of measured data representing a system in first end second states. In the 
embodiment* described herein, fee inversion method estimates the changes of the 
elastic material properties of a region of the Earth containing a hydrocarbon reservoir 
20 due to production or removal of hydrocarbons. The techniques can be based on both 
isotropic and anisotropic models, provided that the reflections may bo expressed 
linearly. The inversion method gives estimates of the changes m the parameters of the 
model, in addition to me corresponding uncertainty. The solution is obtained by 
combining the information provided by difference data with knowledge obtained pnor 
25 to inversion. The solution is therefore more robust and less vulnerable to mstabiht.es. 
Woddng on difference data rather than inverting the measured data prior to taking the 
difference is advantageous with respect to uncertainty estimation in that it allows a 
correct quantitative statistical treatment of the uncertainties, eliminating the need for 
qualitative interpretation. 



15 



30 



To assess the uncertainty of inversion results, the inversion process is cast in a statistical 
setting The solution of an inver se problem is not limited to a single best-fitting set of 
model parameters, but also characterises the uncertainty of the inversion results. A 
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Bayesian setting is chosen for the inversion, although other methods. e.g. least squares, 
can also be used to solve the inversion problem, hi a Bayesian setting it is possible to 
combine available prior knowledge with the information contained in the measured 
data. The solution of a Bayesian inverse problem is represented by foe posterior 
5 distribution, which addresses all questions of nonuniqueness and uncertainty. In 
particular, a Gaussian posterior distribution is completely characterised by a posterior 
expectation and a posterior co-variance. 

A model of foe earth for foe region under consideration is defined, characterised by 
10 parameters describing tte elastic material proparties. The embodiments of the preneiLt 
invention implement a forward modelling operator, the details of which shall now be 

described. 



25 



As described above, an isotropic elastic medium may be fully described using a set of 
three parameters. The embodiments described below adopt foe P-wave velocity o(*. i), 
the S-wave velocity ftx, 0 and foe density p(x, 0 «» parameterisation variables, 
where x is the lateral location and t is the two-way vertical seismic wave travel time. 
An alternative parameterisation which may be used takes the acoustic impedance 
Z„ =(*>, foe shear impedance Z s -jfr. and foe density p as foe parameterisation 
" variables instead. 

The weak contrast reflectivity function r for PP reflections can be written as 

where B is the reflection angle, a a =(l + tau 3 *)/2, =^K/" ^ 
a _ (1 _ 4(£ /«)» sin 2 0)12 . The values of a„ , a p and a„ are defined in 
with a known prior background model of the region under consideration, 
incidence reflections, this reduces to 

KR«.0)=~mZ J >(x,0. 



|0091D25|i09,y3ao^0»iO3Sg»;?t| 



03-JRN-2004 16:37 FROM tlarks&Clerk 

MAC Folio No PS2953GB 



TO PATENT OFFICE 



F.43 



* 



10 



The techniques described herein ate with reference to EP reflections, but the inversion is 
equally applicable to other types of reflection, including PS and SS ructions, and in 
general to all linear expressions of reflectivity. 

A model parameter vector m can be defined as 

m(x,r) «[lna(x,/)uln>?(st»0.biMx>0] r » 
where T denotes transpose such that m(x,f) is a column vector. A further row vector a 

is defined as 

such that the reflectivity function given above may be defined by the dot product 

•a 



set of 
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25 



30 



A discrete representation of the reflectivity r in a given time window for a given 
reflection angles can then be written as a vector 

where A W is a sparse matrix defined by and ut(x) is a discrete representation of 

the model parameters in location x. 

The seismic data a can be represented by the conventional model 

s(x,f,0) * H'»«W V " *• ^ + 4* • 
where w is the seismic wavelet and . is an enor torn, The wavelet may to: angle- 
dependent, out independent of the lateral position x. The wavelet is assumed to be 
nonary within a limited target window. A seismic angle gather at location * can 

therefore be written, using all of the above, as 

s(x) - WA(x)Dm(x) + e(x). 
where w is a matrix representation of the angle-dependent wavelets, A(x) is the sparse 
^ defined by a above, and D is merely a differential operator representing partial 

more compact notation by defining the forward modelling operator G(x> - W to 
give a description of a seismic angle gather at location x of 

s = Gn + e. 
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Seismic data is obtained ftom the region under consideration on two separate occasions, 
shown at step 1 in the flow diagram of Figure 1. Difference data is formed by 
subtracting a first (baseline) seismic dataset si from a subsequent (repeat) seismic 
dataset ffc shown at step 2 of Figure 1 . The method described herein may be applied to 
both the difference of pre-stack seismic data, and the difference of folly stacked or 
partially stacked data. The seismic difference data is discrete and is represented on a 
seismic grid. A seismic difference vector d is defined as a collection of seismic 
difference data from a set of grid locations, and is represented as d = s 2 - si. 

The spatially distributed model parameters are represented on a grid covering the region 
under consideration, using a model parameter vector m containing a collection of the 
model parameters from a set of grid locations, such as those defined above in terms of 
tb* P and S-wave velocities and density. The model grid should covht the region of 
seismic data collection (the seismic grid), but the two do not need to entirely coincide. 
M most cases, however, coinciding grids would be a preferable choice. ^ The 
embodiments described herein refer to coinciding grids. The step e" 
parnmeterised model is shown at step 3 of Figure 1 . 



defining 



20 A model parameter change vector is defined as the difference between the model 
parameter vector m, corresponding to the repeat dataset and me model parameter vector 
m, corresponding to the first dataset. that is o = -* - m, , corresponding to step 4 of 
Figure 1. 



25 



30 



The model parameter change vector 6 may be defined (tor constant - 
parameterisatian above, as 

where the indexes 1 and 2 refer to the first and second seismic data gathers respectively. 
For zero-incidence reflections, this reduces to 

where Z F » « a* is the acoustic impedance at time /. 
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A linear relationship between the model parameter change vector and the seismic data 
difference vector can be written in matrix-vector notation as 

d «» GS + e, 

where G is the forward modelling operator G(x) defined above, and e is the error term 
being related directly to the difference data d. The formulation of this inversion 
problem represents a new approach to the inversion of time-lapse seismic data. 

The error term e is predominantly a consequence of seismic noise which may be 
chacacterked with an enor expectation vector and an error covariance matrix 1^ The 

•i-r« £_ — wBAiv-tv . 







■ 
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correlations between the different elements therein. The error covariance inatn* can be 
e5timated fton, the seismic difference data in regions not influenced by the production 
under consideration, which is not possible unless data are differenced. The error 
expectation and error covariance are determined prior to the inversion. 

The knowledge and uncertainty about the model parameter change vector 8 prior to the 
inversion are characterized via a model parameter change expectation vector M and a 
model parameter change covariance matrix £3. The model parameter change 
expectation vector and the model parameter change covariance matrix axe specified 
prior to inversion, and can be determined by analysis of the effects of fluid substaut.cn 
^d pressure changes due to production. In regions not affected by production, the 
model parameter change expectation is «*o (no change is expected). The model change 
parameter covariance matrix specifies the variance for each of the elements fa**/*). 
ln(pVP0, mWPt) * ^ model parameter change vector, and the correlations 
therebetween. 

Explicit analytical expressions are calculated for a solution in which the prior model is 
combined with the information provided by foe seismic difference data by forming a 

posterior distribution, shown at step 5 of Figure 1. The solution is represented via an 
updated model parameter change expectation vector, step 6 of Figure 1, and an updated 
model parameter change covariance matrix, step 7 of Figure 1. The updated model 
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parameter change expectation vector provides an estimator of the change in material 
properties, and thus allows quantitative inferences on the nature of the change to be 
made. The uncertainty is represented via this updated covariance matrix. The analytical 
explicit expressions for the solution provide a computationally fast inversion method 
which allows for the assessment of the likelihood of changes in the physical system, 
shown at step 8 of Figure 1. 




10 



15 



20 



The statistical properties of the model parameters are assumed to be Gaussian, although 
any other statistical distribution, eg. the Cauchy distribution, may be used. If each 
single parameter is assumed to be Gaussian defined by an expectation value and a 
variance, then the model parameter change vector 6 is Gaussian, defined by the model 
parameter change expectation vector and the model parameter change covariance 
matrix: 

where n m is the dimension of the model parameter change vector, and |i« end S* ate as 
defined above. 

The statistical properties of the noise are also assumed to be Gaussian. If each single 
seismic noise sample is assumed to be Gaussian denned by an expectation value and a 
variance, then a vector of noise samples is Gaussian defined by an error expectation 
vector and an error covariance matrix as described above: 

dlS-JV./Go.E,), 







1 


El 



25 described above. 

These two distribution models therefore imply mat the distribution for the seismic 
difference data is also Gaussian; 

30 with a difference expectation vector m - Go, and difierence covariance matrix 
S,=GE s G r +L„. 



l0O91oa:j09.?=aaife.04:.03i3»:ll 
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Suppressing the location parameter x, the joiat distribution for the model parameter 
change vector 5 oand the seismic difference data vector d is then: 



'd J 



lie noise. 



Under the described GausBian assumptions for the prior model and the seii 
the solution is represented via die multi-Gaussian distribution: 

defined by the new combined expectation and covariance. In this case, the expectation 
vector, represents an .optimal solution for jhe. 4D inversion problem. In Bayesian 
terminology, this is the posterior distribution wim a posterior expectation of 

and a posterior covariance of 



The posterior expectation is an optimal estimator for o calculated from the 
data d, such diet 

while the uncertainty is represented by the posterior covariance matrix. 



20 
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The posterior estimator Sis an estimator for btfaa/ai), ln(pVP0 and Info/pi), hi the 
poststack case, given aero incidence reflections, $ is an estimator for ta(ZpyZp,i). 
Therefore 5 gives a quantitative estimate of the change in material properties, and the 
corresponding posterior covariance matrix gives its uncertainty. These two quantities 
may then be used to infer changes in production in the reservoir, and other related 
reservoir properties. 

The method of the above embodiment will now be illustrated byway of example. 

A reservoir model ia defined by acoustic impedance models for the baseline and repeat 
seismic surveys, illustrated in Figures 2a and 2b. The reservoir is originally saturated 
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by oil t but during production, regions of the reservoir arc water-flushed. The drainage is 
modelled by m increase of the acoustic impedance. The percentage change in acoustic 
impedance is illustrated in Figure 2c. where the brighter-coloured region indicates the 
region where the physical properties have changed. 

i 

For each of the baseline and Tepeat models, a seismic forward modelling is performed 
by convolution. The wavelet is a. Ricker wavelet with a 25Hz centre frequency and 
normalised amplitude. The modelled results 81 and s 3 are shown in Figures 3a and 3b 
respectively. A combination of coloured noise ei and white noise ej of the form 
e=s, *e, +e 9 was added to the difference data d. The coloured noise represents 
coherent "noisi, for example a non-repeatability of the source energy between the 
baseline and repeat surveys. The white noise is Gaussian with a variance of <r\ . The 
seismic difference was simulated with signal-to-noise ratios ofS/N=10 5 and 2, shown in 
Figures 3c and 





To test uie inversion method, a stationary prior distribution is assumed. i.e. 




I- 



which coiresponds to_a constant acoustic impedance between surveys. A atationary 
covariance function with a* ■= 0.039 is also assumed, corresponding to a 0.95 prior 
model interval of ±8% acoustic impedance change, which widely includes me modelled 
increase of about 2-4%. The noise covariance and the wavelet used in me inversion are 
consistent with the model. 

The result of invention of the difference data with a S/N ratio of 10 s is shown in Figures 
4a to 4c. where Figure 4a is the true model, Figure 4b is a waveform representation of 
the seismic difference data, and Figure 4c shows the posterior mean solution. With this 
low noise level the acoustic impedance change is retrieved almost exactly, with low 
uncertainty. 



30 The result of inversion of the difference data with a S/N ratio of 2 is shown in Figures 
5a to 5c. where Figure 5a is the true model. Figure 5b is a waveform representation of 



09-JRN-2004 16! 33 FROM Marks&Clerk 

M&CPbUoNoH2353GB 



TO POTENT OFFICE 



P. 49 




10 



15 



20 



25 



13 

the seismic difference data, and Figure 5c shows the posterior moan solution- The 
solution is smoother and has higher uncertainty that that of Figure 4c. 

The uncertainties of the inversion results at CDP 115 are shown in Figure 6a fbr the low 
noise case (S/N-10 5 ), and in Figure 6b fbr the Wghnoise case (S/N-2). The confidence 
region in Figure 6b is much wider than in Figure 6a, as would be expected. 

Based on the uncertainty bounds of tho inversion results, it is possible to provide 
quantitative probabilities of different reservoir states, e.g. drained, unchained, amount of 
vertical change of a hydrocarbon contact. The probabilities can be used in risk value 
estimation and commercial decisions. In this example, Figure 7 illustrates a probability 
density map of drained regions (witb S/NH2), where darker colouring indicates a higher 
probability of having been drained. This corresponds very well to the modelled case of 
Figure 2c. 
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The data inversion methods described above may be embodied in a program for 
controlling a computer to perform the inversion. The program may be stored on a 
storage medium, for example hard or floppy discs, CD or DVD-recordable media or 
flash memory storage products. The program may also be transmitted across a 
computer network, for example the Internet or a group of computers connected together 

» • 

in a LAN. 

* 

The schematic diagram of Figure 8 Ulustrates a central processing unit (CPU) 13 
connected to a read-only memory (ROM) 10 and a random access memory (RAM) 12. 
The CPU is provided with measured data 14 and model parameters 16 via an 
input/output mechanism IS. The CPU then performs the inversion on the provided data 
in accordance with the instructions provided by the program storage (1 1) (which may be 
a part of the ROM 10) and provides the output, i.e. the updated model parameters and 
uncertainties 17, Via the input/output mechanism 15. The program itself, or any of the 
inputs and/or outputs to the system may be provided or transmitted to/from a 
communications network 1 8, which may be, for example, the Internet. 
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CLAIMS: 



1. a method of processing seismic data representing a physical system, the method 
comprising: 

forming! a difference between first and second seismic data representing the 
system in first and second states, respectively; and 

inverting the difference in accordance with a parameterised model of the 
physical system to obtain changes in the parameters of the modeL 



10 2. A method as claimed in claim l f wherein the first and second states 
physical system at different times* 



15 



3. 



in a 



A method as claimed in any preceding claim, wherein the inversion is performed 



4, 



A method as claimed in claim 3, wherein the inversion is a least squares 



inversion. 



20 



5. 



A method as claimed hi claim 3, wherein the inversion is a Bayesian inversion. 



6. A method as claimed in claim 4 or 5, wherein a first part of the solution to the 
inversion is a posterior expectation. 

7. A method as claimed in claim 6> wherein the posterior expectation corresponds 
25 to changes in the parameters of the model. 

8. A method as claimed in any one of claims 4 to 7, wherein a second part of the 
solution to the inversion is a posterior covariance. 



30 9. A method as claimed in claim 8. wherein the posterior covariance corresponds to 
the uncertainty in the posterior expectation. 



1 po9io25:.:ro.9>:J mm^^m^i 
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10. A method as claimed in any preceding claim, wherein changes in tho parameters 
of the model follow Gaussian statistics. 



11. A method as claimed in any preceding claim, wherein the difference has an 
5 associated error model with a corresponding error term. 



12. A method as claimed in claim 11, wherein the error term corresponds to the 
difference between (i) the difference between the first and second measured data, and 
(il) the changes in the parameters of the model when operated upon by a forward 

10 modelling operator. 

13 . A method as claimed in claim 1 1 or 12, wherein the error term e is related to the 
difference d between the first and second measured data and the changes 6 in the 
parameters of the model when operated upon by a forward modelling operator G via the 

15 matrix-vector expression d = <S8 + «- 

14 . A method as claimed in any one of claims 11 to 13. wherein the statistical 
properties of the error term is described via an expectation and a oovariance. 



20 15. 



25 



A method as claimed in claim 14 wherein the error term fellows Gaussian 



30 



statistics. 

16. A method as claimed in any one of claims 1 1 to 15 when dependent on claim 6, 
wherein the cm* term is used to determine the posterior expectation. 

17. A method as claimed in any one of claims 1 1 to 16 when dependent on claim 8, 
wherein the error term is used to determine the posterior covariance. 

* 

18. A method as claimed in any one of claims 11 to 17, wherein the physical system 
has a part which is the same in the first and second states, measurements of this part 
being used to determine statistical properties of the error term. 
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19. A method as claimed in any preceding claim, wherein the prior knowledge about 
the changes of the parameters of the model is defined in a prior model. 

20. A method as claimed in claim 19, where the prior model is defined via a prior 
expectation and a prior covariance. 

21. A method as claimed in claim 20, where the prior model follows Gaussian 




statistics. 



of claims 19 to 21 when dependent on claim 1<S, 



10 22. A method as claimed in any one 

where the prior model for the changes of the parameters of the model is used to 

determine the posterior expectation. 

23 A method as claimed in any one of claims 19 to 22 when dependent on claim 17, 
15 where the prior model for the changes of the parameters of the model is used to 
determine the posterior co variance 



24. A method as claimed in 
reflection data. 



any preceding claim, wherein foe seismic data are 



20 



25 A method as claimed in any- preceding claim, wherein the parameters of the 
model are any complete set of elastic properties that can be used to formulate a linear 
expression of reflectivity. 



4 r - 



claim 25, wherein the parameters of the model are the P- 



25 26. A method as claimed m 

wave velocity, the S-wave velocity, and the density. 

27. A method as claimed in claim 25, wherein the parameters of foe model are the 
acoustic impedance, foe shear impedance, and tho density. 

30 , 

28. A method as claimed in any preceding claim, wherein foe physical system 

includes a region containing a hydrocarbon reservoir. 
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29. A method && claimed in any preceding claim, wherein the obtained changes in 
the parameters of the model are used to assess where changes in die physical system 
have occurred between the first and second states. 

5 30. A method* as claimed in claim 29 when dependent on claim 28, wherein the 
obtained changes in the parameters of the model ate used to assess the change in 
hydrocarbon content in the hydrocarbon reservoir. 

31, A method as claimed in any preceding claim, wherein the obtained changes in 
10 the parameters of the model are used to assess the probabilities of changes in the 
physical system between the first and second states. 



physical system are used to generate a probability density map of changes in the 
IS physical system. 

* 

33. A program for controlling a computer to perform a method as claimed in say 
one of claims 1 to 32. 

20 34. A program as claimed in claim 33 stored on a storage medium. 



35. Transmission of a program as 

network. 



in claim 33 across a communications 



25 36. A computer programmed to perform a method as claimed to any one of claims I 
to 32. 



37. An apparatus for processing data representing a physical system, the apparatus 
comprising: 

30 means for forming a diffe*nce between first and second measured data 

representing the system in first and second states, respectively; and 

means for inverting foe difference in accordance with a paiameterised model of 
the physical system to obtain changes in the parameters of the model. 



I* 
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ABSTRACT 

Processing seismic data representing a physical system 



10 



A method is provided for processing seismic data representing a physical system, an 
as seismic reflection data representing a region of the Earth. The method forms a 
difference (2} between first and second measured data (I) representing the system in 
fcst and second states, respectively, for example at different times. The difference is 
inverted (5) in accordance with a parameterised model (3. 4) of the physical system to 
obtain changes (6. 7) in the parameters of the model. 



Figure 1 
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